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A THERMODYNAMIC CHART FOR THE MARS ATMOSPHERE 


By David E. Pitts 
Lyndon B. Johnson Space Center 


SUMMARY 


A thermodynamic chart known as a tephigram (temperature /entropy gram) was 
drawn for the Mars atmosphere by using an X-Y plotter and a computer. The teph- 
igram has temperature as the abscissa and the log^ of potential temperature as 

the ordinate. Lines of pressure, temperature, potential temperature , saturation- 
mixing ratio, and equivalent potential temperature are presented so that water- 
vapor-condensation processes may be studied for a range of temperatures from 
80 to 320 K, for a range of potential temperatures from 100 to 390 K, and for a 
range of pressures from 0.1 to 8 millibars. 

The region in which carbon dioxide condensation can occur is shown at tem- 
peratures of 150 K or less and is represented by an area on the chart. A min- 
imum of assumptions was made so that maximum use could be made of the chart. 

These assumptions , that the surface pressure is 3 millibars and the atmosphere 
is 100 percent carbon dioxide, are consistent with analyses of the Mariner IV 

-3 

data. Trace amounts of water vapor (approximately 1 * 10 centimeter precipi- 
table) were assumed. 


INTRODUCTION 


The fact that so many disciplines are involved in the interpretation of 
phenomena in the Mars atmosphere increases the difficulty of coordinating a 
self-consistent analysis of the state of the atmosphere. Thus, the purpose of 
this report is to present a meteorological chart for the atmosphere of Mars that 
will aid investigators in the interpretation of in situ measurements (e.g.. 
Viking), remote measurements from orbiters or from flyby missions, astronomical 
observations, and meteorological theories. 
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SYMBOLS 




specific heat at constant pressure for carbon dioxide gas, calories/ 
mole-kelvin 


c^ constant 

dq change of heat in a unit mass of gas 
e . saturation vapor pressure for water ice, millibars 

S X 

e saturation vapor pressure for water, millibars 
sw 

f factor to account for nonideology of gas 

2 

G acceleration of gravity for Mars, 375 centimeters/second 



L latent heat of vaporization for water, calories/gram 
w 

2 

M col umn ar mass, grams /centimeter 

m molecular weight of ambient atmosphere 

m, molecular weight of carbon dioxide gas, HU. Oil grams /mole 
d 

m molecular weight of water vapor, 18.016 grams /mole 
w 

n number of approximations 

P pressure , millibars 

P surface pressure, millibars 

o 

P pressure at altitude z, millibars 

z 

P^ same pressure lower than Pg, millibars 
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some pressure higher than P^, millibars 

Q = 8(1/T)" 1/K 

R, universal gas constant, 1 . 98583 calories /mole-kelvin 
a 

r mixing ratio, grams /kilogram 

r average mixing ratio from P to P , grams /kilogram 

o z 

r saturation-mixing ratio, grams/kilogram 

S 

T temperature , kelvin 

T average temperature between P^ and P^, kelvin 

Tj dewpoint temperature, kelvin 

1 3 

o specific volume — , centimeters /gram 

3 constant from reference 9 used in equation (30) 

AZ geopotential altitude between P^ and P^, meters 

n constant from reference 9 used in equation (30) 

0 potential temperature based on surface pressure of 8 millibars and an 

atmosphere of approximately 100 percent carbon dioxide, kelvin 

6g equivalent potential temperature based on a change of phase of water in a 
100-percent-carbon-dioxide atmosphere, kelvin 

£ constant from reference 9 used in equation (30) 

p atmospheric gas density, grams /centimeter^ 

4> entropy per unit mass 

CONSTRUCTION OF THE MARS TEPHIGRAM 

Saucier and Elliott (ref. l) developed a thermodynamic chart or tephigram 
for the Earth atmosphere that has proved extremely useful. A tephigram for the 
Mars atmosphere has been developed that will aid investigators in interpreting 

3 



the maximum amounts of water vapor from the range of surface temperatures ob- 
served by the Mariner IX infrared radiometer (ref. 2) of 195 K (morning termina- 
tor) to 250 K (local noon) and in interpreting temperature structure deduced 
from S-band occultations (ref. 3) and from temperature sounders like the infra- 
red interferometer spectrometer (ref. 4) or from entry probes such as the Viking 
lander. A copy of the Mars tephigram is included at the end of this report. 

Most of the major features of the tephigram are common to other thermody- 
namic charts; however, the tephigram was chosen because the right angle between 
the isotherms (temperature T) and dry adiabats (potential temperature 
6 lines) makes it easier to detect variations in the temperature lapse rate. 
Area is also proportional to energy, as is true of any true thermodynamic chart. 

Assumptions consistent with Mariner IV and Mariner 1Y data were used (i.e., 
a surface pressure of 8 millibars and a 100-percent-carbon-dioxide (CO^) atmos- 
phere, with trace amounts of water ( H^O ) ) . A. J. JCLiore et al. (ref. 3) found 

a range of surface pressures between 2.9 and 8.3 millibars and elevation differ- 
ences as great as 13 kilometers from Mariner IX occultations at S-band frequen- 
cies. Recent models have used 50 to 100 percent CO^ corresponding to surface 

pressures of 15 and 5 millibars, respectively (ref. 5) • However, with lower 
surface pressures being corroborated more recently, there is increasing use of 
100 percent C0^ (ref. 4) primarily because no other gas has been detected in 

large quantities (ref. 5)* 

Lines of pressure, temperature , potential temperature, saturation-mixing 
ratio, and equivalent potential temperature are presented so that water-vapor- 
condensation processes may be studied for a range of pressures from 0.1 to 
8 millibars, for a range of temperatures from 80 to 320 K, for a range of po- 
tential temperatures from 100 to 390 K, for a range of equivalent potential 
temperatures from 200 to 400 K, and for a range of saturation-mixing ratios 

r (0/00; i.e., grams per kilogram) from 10~^ 0/00 to 10^ 0/00. 
s 

A line representing the equivalence of the saturation vapor pressure of 
CC^ and total, pressure is drawn on the tephigram to indicate the area where CO^ 

condensation can occur. A temperature of 150 K or less is necessary at the 
surface; at 0.2 millibar, 126 K or less is required. This temperature is a 
necessary but not sufficient condition for condensation because small droplets 
require higher vapor pressure for equilibrium than does the flat surface as- 
sumed in the vapor pressure equations. Because CO^ is the major constituent 

of the atmosphere, condensation of C0^ will not effectively change the amount 

of C0p (i.e., mixing ratio); however, it will tend to stabilize the temperature 

at the freezing point of CO g and cause a reduction in surface pressure. 

Using a computer program, approximately 20 000 data points were generated 
for the tephigram. These data were processed by an X-Y plotter that draws lines 
to an accuracy of five significant decimal digits. The variability with temper- 
ature of specific beat at constant pressure and the latent heat of vaporization 
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should be accounted for with such accurate plotting and computer capability; 
this was accomplished with second-order curve fits on temperature. All quanti- 
ties were plotted relative to the ordinate log^g ® and abscissa T; that is, 

no lines were plotted relative to the pressure lines or the saturation-mi xing- 
ratio lines. 


GRID STRUCTURE 


Isotherms are the vertical yellow lines entered for each 1 K; the heavy 
lines are for each 10 K in temperature. The isotherms are labeled at the lower 
ends. 


Dry adiabats or isopleths of potential temperature and of entropy are the 
horizontal yellow lines that are drawn in the same scheme as the isotherms and 
labeled in kelvin. The vertical scale is log^g 9, which is directly propor- 
tional to log e 9 where 9 is potential temperature in kelvin. Thus, because 
specific entropy (entropy per unit mass) 


= C pd lo8 e 6 * C 1 


(1) 


where C ^ is specific heat at constant pressure for CO^ gas and is a 

constant, and "because the horizontal scale is proportional to temperature, the 
chart is a temperature /entropy diagram or tephigram. 


Area on the tephigram (ref. 1) shows that one form of the first law of 
thermodynamics is 


dq 



m 


a dP 


(2) 


where dq is change of heat in a unit mass of gas, m is molecular weight of 
ambient atmosphere, P is pressure, and a is specific volume 1/p. B cause 


a - 


V 

Pm 


(3) 
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where R^ is the universal g as constant, equation (2) may be rewritten as 


C R 

IP* = d(ln T) - — d(ln 
T m m 


P) 


(U) 


Both terms on the right-hand side of this equation are exact differentials , so 


/j 1 -/ 4 ** 0 


(5) 


Potential temperature is defined as the temperature corresponding to the 
surface pressure P q . By the integration of equation (U), from P Q to P 

(6 to T), using equation (5) , one obtains 


e = 



( 6 ) 


where P « 8 millibars for Mars. Then 
o 


R j Rj 

In 0 * In T + 7T“ In P - In P 
C . o C , 
pd pd 


(7) 


Taking the differential gives 


d(ln 0) = d(ln T) - — d(ln P} (8) 

Pd 

Substituting equation (8) into equation (U) gives 


dl 

T 


C . d(ln 9) 
Pd 


( 9 ) 
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Thus, the work done dq in a cyclic process is 



'pd 


/. d(ln 6) 


do) 


Therefore, because T is the abscissa and log 0 is th r lina. area on the 
tephigram represents energy. 

Isobars are the green curves sloping upward to the rigl.t, labeled in milli- 
bars, and drawn in 0.1-millibar increments from 8 millibars to * millibar and 
in 0.01-millibar increments from 1 millibar to 0.1 millibar. Tney were calcu- 
lated by using the following equations. Equation (6) in log form is 


so 


In 0 = In 



( 11 ) 


log 


10 


0 = 0.U3U29UU8 In T 



( 12 ) 


Isopleths of saturation-mixing ratio r g are the green lines that slope 

upward to the left and are not greatly curved. Saucier and Elliott (ref. 1) 

define r ar 
s 


r 

s 


m f e 
w w sw 


m 


,/p - f e ~\ 

l\ w sw J 


(13) 


where m w is the molecular weight of water vapor, f is a factor to account 
for nonideology of gas, e is the saturation vapor pressure for water, and 

FW 

ittj is the molecular weight of CO^ gas. For the low pressures involved here. 
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f = 1 may be assumed and still be accurate within ±0.02 percent (ref. 6). 

W 

Thus, equation (13) may be simplified to 


r 

s 


m 

w 




(lU) 


However, for plotting purposes, r had to be a function of T and log 0. 

S XU 

Equation (6) may be solved for P 



and then substituted into equation (lU) giving 


m e 
V sv 


m p * 

s 

-Ed 



m d 

Kf) -vj 


Solving for 8 gives 


8 


T 



m.r 8 
d s 

+ m.r \e 
d s J sw 



(15) 


(16) 


(17) 


or 


lo8 10 


8 » O.U3U29UU8 


In T + 


'pd 



(18) 
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The saturation vapor pressure e gv that was used in equations (l6) to (16) was 
taken from Saucier (ref. 7). For water 


e 

sw 


= 6.11 x 10 [7.5T/(23T.3*T)] 


For ice 


e . 
01 


6.11 x 1 c [9.5T/(26 5 .0 + T)] 


(19) 


( 20 ) 


where e ^ is saturation vapor pressure for water ice and T i in degrees 
Celsius . 

Isopleths of equivalent potential temperature, which are approximate 
ps eudoadi abate , are the curved green lines that are concave toward low tempera- 
ture and potential temperature and are entered for each 10 K from 200 to U00 K. 

The values of equivalent potential temperat ure 0- are found by following the 

h 

pseudoadiabat to the left border of the chart and reading the value m the po- 
tential temperature scale; that is, 0 approaches 0 at low temperature. 

The 0 lines were drawn by using the expansion of the Rossby equation (ref. l) 



( 21 ) 


where is the latent heat of vaporization for water, and solving for 

6 = f(T). Since r = f(0,T) in equation (l6), 0 must be found by taking 

s 

successive approximations of 0 in the form 


f(e) * 0 (l 


(■ 



( 22 ) 


until f(0) « 0. Successive approximations are found by Newton's method 
(ref. 8). 
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(23) 


e 


n+1 



for n approximations . Five such approximations were sufficient for each data 
point when solutions for the 6g line were started at low temperature and 

worked to higher temperature 1 K at a time. The preceding solution is used each 
time as the first guess to the next point. 

Taking the derivative of equation (23) with respect to 9 gives 


f'(e) = 



6 srH 1 ♦ <*, 




(2U) 


where 


C 



(25) 


K * 



(26) 


The derivative of r with respect to 0 in equation (lb) is 

s 


where 



(27) 



( 28 ) 
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During the calculations of the P, r , and 6_ lines* C . was assumed 

* s E pa 

to vary with temperature according to the formula 


c = 6.6367 ♦ 1.396 * lO'-^T ♦ 2.0U15 x io~ 5 t 2 (29) 

pd 

which was obtained by performing a curve fit of the form 


C . = n ♦ BT ♦ CT 2 (30) 

pd 


on data from McBride et al. (ref. 9)» where n» 8* and £ are constants. The 
value of L w was assumed to vary in a manner similar to that in equation (30). 

In this case, the curve fit was performed using data from List (ref. 6). 


L w = 816.9^32937 - 1.00 509 12 37T + 0.0007°6001T £ 


(31) 


Performance of the thickness evaluation was similar to that done by 
Saucier and Elliott (ref. l) with modifications for the Mars atmosphere. 


AZ = 50.3777T log 


a 


(32) 


where the depth of the layer AZ in geopotential meters is proportional to the 
log e of the pressures at the lower P and upper P^ boundaries and the mean 

temperature of the layer T. 

Hie carbon dioxide equilibrium saturation line is the very heavy line on 
the left-hand side of the chart. It was drawn by assuming that the pressure of 
the atmosphere (100 percent C0„ ) was equal to the saturation vapor pressure for 
C0 2 (ref. 10). Z 

(, . 9 . 9oa2 ) 

P « 1.333225685 * 10' (33) 


where T is in <* grees Celsius. The region to the left of this line represents 
the saturation r gion for a plane surface of solid CO 2 » This statement, of 
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course* does not mean that saturation must occur in this region. For example, 
when small particles of CO^ are involved, saturation will not occur until a 

lower temperature is reached. When C0_ sublimes, 137.037 cal/g are involved 
(ref. 11). d 


GRAPHICAL OPERATIONS 


The processes briefly described in the following sections are some of the 
more common operations performed with the tephigram. A more detailed treatment 
can be found in reference 7. 


Adiabatic Processes 

Air heated near the surface ascends adiabatically (6 is conserved); that 
is, the pressure and temperature can be found by following the 9 line cor- 
responding to the original pressure and temperature . If the surface pressure 
is 8 millibars and the surface (air) temperature is 230 K, rising air stays 
at a constant 9 = c30 K. Thus, when the pressure is U.S nillibars , the tem- 
perature is 197 K. 


Saturation 

The amount of water necessary for condensation at any pressure and temper- 
ature is given by the corresponding r g as read from the chart. For example, 

for a pressure of U.5 millibars and a temperature of 197 K, the saturation- 
mixing ratio r g is 0.1 0/00. If the mixing ratio r = r g , saturation occurs 

and the parcel no longer ascends adiabatically. Instead, the parcel follows a 
line corresponding to the pressure and temperature where r = r (the base 

ii* S 

of the cloud). If the surface pressure is 8 millibars, the surface temperature 
is 230 K, and the mixing ratio is 0.1 0/00, the parcel is lifted adiabatically 
and conserves moisture r at U.5 millibars pressure and 197 K. Above this 
altitude, pressure and temperature values can be read from the 0g line of 
232 K because the processes are such that 0^ is conserved. 


Isobaric Processes 

Isobaric cooling from the original pressure and temperature with P re- 
maining constant without addition or removal of moisture (r is constant) re- 
sults in saturation when r * r g . This process defines the dewpoint temperature 

T.. Isobaric cooling results in stability or inversion and will produce fog if 
d 

the moisture content r i3 sufficiently high. If the surface pressure is 
8 millibars, the surface temperature is 230 K, and the water vapor is 0.1 0/00; 
then the dewpoint temperature T^ is 202 K. Isobaric warming near the surface 
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causes instability. This condition produces vertical motion and is commonly 
thought to produce clouds on Mars. 


Stability 


Unsaturated air is stable or unstable according to the relationship of the 

30 30 

temperature lapse rate to that of the 9 lines : - 0 is neutral , rr- < 0 

30 *** 

is stable, and ^ > 0 is unstable. Saturated air is stable or unstable ac- 
cording to the relationship of the temperature lapse rate to that of the 0^, 


lines. Thus, it is possible that an air column could be stable for unsaturation 
and unstable upon saturation. This particular situation is called conditional 
instability. Thus, above saturation, the parcel of air would be in free con- 


vection (at the '■»vel of free convection): 


unstable, and 




is stable. 


i£ 

3P 


30 e 

— r— is neutral. 




3P 


is 


Effects of Diurnal Heating 

As morning comes, the nighttime inversion is removed and the ground tem- 
perature becomes progressively higher as heating continues. The air column is 
modified to an adiabatic state in the lower regions while the higher regions 
remain undisturbed. The area between the final adiabatic lapse rate, the sur- 
face pressure curve , and the original sounding represents the energy added by 
insulation. If the heating from below is intense enough, condensation may oc- 
cur. The point at which condensation occurs is the convective condensation 
level. 


Mighttime Fog Due to Inversions 

When rapid cooling near the surface occurs, J ' mixing ratio may be reached 
by the saturation-mixing ratio in the resulting pi- .sure temperature sounding. 
This condition results in nighttime fog. 


Columnar Mass of Water 


The amount of H^O in a column 


(grams per square centimeter) is equal to 


M = 1000 



r(0/00) p(z)dz 


(3U) 
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where M is the columnar mass of H^O, p is the atmospheric gas density, and 

r is in grams per kilogram. However, if the average mixing ratio r from the 
surface pressure P to the pressure P at altitude z is substituted for 
r, then z 



(35) 


where G is the acceleration of gravity for Mars. To find the amount of pre- 
cipitable H^O in centimeters, divide by the density of water (l g/cm ). 

CONCLUDING REMARKS 


The tephigram of the Mars atmosphere presents a means of performing de- 
tailed graphical operatic ns such as adiabatic processes, saturation, isobaric 
processes (such as radiative cooling), stability changes with condensation, and 
other processes without the aid of costly computer analysis. Thus, analyses of 
the Mars atmosphere may be facilitated by the Mars tephigram as studies of the 
Earth atmosphere have been by the tephigram of Saucier and Elliott. 


Lyndon B. Johnson Space Center 

National Aeronautics and Space Administration 
Houston, Texas, March 6, 197 
951-16-00-00-72 
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